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£h Abstract 
On 

It is known that a higher concentration of free carriers leads to a higher oxide growth rate in the 
thermal oxidation of silicon. However, the role of electrons and holes in oxidation chemistry is not 

"o 

clear. Here, we report real-time second-harmonic-generation data on the oxidation of H-terminated 
(111) Si that reveal that high concentrations of electrons increase the chemical reactivity of the 
outer-layer Si-Si back bonds relative to the Si-H up bonds. However, the thicknesses of the natural 
oxides of all samples stabilize near 1 nm at room temperature, regardless of the chemical kinetics 
i-^j of the different bonds. 
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The passivation of Si substrates and the chemical stability of passivated surfaces are both 
significant for Si-based device technology. [1J The most common way to passivate silicon 
surfaces is to terminate them with H by dipping them into HF or NH 4 F solutions. This 
process strips the oxide layer and caps the outermost orbitals of the Si substrate with H. 

The chemical stability of these surfaces depends on the reactivity of two types of bonds: 
the outermost Si-H bonds and the Si-Si back bonds to the underlying layer. Determining the 
individual chemical reactivities of these bonds is essential for understanding and controlling 
chemical processes on H-terminated silicon surfaces. Under ambient conditions the domi- 
nant chemical process is oxidation. The oxidation rate depends on environmental conditions 
such as temperature, pressure, humidity, and the possible presence of reacting agents as 
well as intrinsic sample properties such as surface morphology and type and density of free 
carriers. [2HU The effect of carrier concentration on oxidation at elevated temperatures has 
been characterized by spectroscopic ellisometry (SE), [SJ E] surface differential reflectance 
spectroscopy (SDR) [7], Auger electron spectroscopy (AES)[8J, and X-ray photoelectron spec- 
troscopy (XPS)[Qj, all of which probe oxide thickness. The oxidation rate has been found to 
increase with increasing carrier concentration. 

The studies cited above determine the macroscopic behavior of oxidation, but do not 
provide information about the reactivities of the individual types of bonds. Here, we report 
results obtained by second-harmonic-generation anisotropy (SHGA) and SE measurements 
that address the chemical reactivities of the different types of bonds that are present in H- 
terminated (111) Si. The bond sensitivity of SHGA is apparent when the data are analyzed 
using the anisotropic bond model (ABM) (TUJ [TT] of nonlinear optics. This approach shows 
that the radiated SHG signal is maximized when a particular bond is aligned parallel to the 
driving field, i.e., when the acceleration of that bond charge is itself maximized. As samples 
are rotated, the back bonds realize this condition sequentially, allowing the respective hy- 
perpolarizabilities to be extracted from SHGA data. [12] Therefore, changes in the chemical 
nature of these bond are easily detected by SHGA. The SHGA response of the "up" bonds is 
isotropic, allowing their chemical reactivity to be inferred from the signal offset. Our results 
show that under ambient laboratory conditions the reactivity of the first-layer Si-Si bonds 
is significantly affected by carrier concentration. Nevertheless, despite these differences in 
early oxidation rates, overall oxide thicknesses evolve similarly and stabilize near 1.0 nm for 
all doping levels. 
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The data reported here were obtained on 10 12 , 10 14 , and 10 16 cm 3 P-doped (lll)Si wafers 
and one 10 18 cm -3 As-doped (111) Si wafer, all purchased from Virginia Semiconductor. The 
samples were cleaned by consecutive 10-min immersions in 80 °C NaOH/H 2 02/H 2 (1:1:5) 
and 80 C HCI/H2O2/H2O (1:1:5) solutions. Native oxides were then stripped and the samples 
capped with H by a 20-min immersion in 40% NH4F. To prevent pitting, the NH4F solution 
was deoxygenated prior to immersion. [T3] Measurements began approximately 90 s after 
NH 4 F immersion, after the surfaces were dried with high-purity N 2 . 

SHGA data were obtained in a reflection geometry using p-polarized light at an angle 
of incidence of 45° for both illumination and detection. The exciting beam was generated 
by a Ti-sapphire oscillator. It consisted of 100 fs pulses centered at 806 nm, arriving at a 
repetition rate of 70 MHz. The SHG response was detected by a photomultiplier every 1° 
as the sample was rotated by 360°. Further details are provided elsewhere. |12j 

Figures la, b, and c show the data for the P-doped samples. All are normalized to the 
level where the SHG responses no longer evolve with time. Where present, the three major 
peaks correspond to the near alignment of one of the three back bonds to the direction 
of the exciting field. As expected, the overall effect of oxidation is to increase the magni- 
tude of the SHG signal at these azimuths, a consequence of the higher electronegativity of 
O relative to H. The evolution clearly depends significantly on carrier concentration. For 
the lowest two carrier concentration the SHG intensity initially increases then decreases to 
its terminal value, whereas for the most heavily doped sample the SHG intensity increases 
monotonically. This striking dependence on carrier concentration is even more evident in 
Fig. Id, which shows the average evolution of the signal at these three azimuths for all four 
samples, including that of the As-doped sample for which n =10 18 cm~ 3 . To relate these 
differences to the chemical kinetics of the bonds, we determined the respective hyperpolar- 
izabilities by fitting the data to the predictions of the ABM using expressions given in [10]. 
The hyperpolarizability results are summarized in Fig. 2. 

The effect of chemistry on hyperpolarizabilities can be understood as follows. Generation 
of SHG by a bond requires that the bond be asymmetric. While an Si-H bond is obviously 
asymmetric, it is less obvious that terminating an up bond with H generates contributions 
from the Si-Si back bonds as well. The reason is chemical induction. The greater electroneg- 
ativity of H removes charge from the Si atom to which it is attached, making the Si-Si back 
bonds asymmetric as well. Thus the SHG contribution of the back bonds tends to track 
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that of the up bond. If the top H is now replaced by OH, the greater electronegativity of O 
results in greater asymmetries and hence stronger signals for both up and back bonds. This 
type of response is clearly evident in Fig. 2 for the lightly doped sample. 

With this background, the striking difference among the samples of different carrier con- 
centrations can now be explained. In all cases, oxidation begins with the replacement of 
H with OH, as described above. The next step is the insertion of the O into a back bond, 
leaving the surface still capped with H. For lightly doped samples this step is slow. But for 
heavily n-type samples the reactivity of the back bonds is much higher, and replacement 
and insertion occur at similar rates. Thus the sequential increase and decrease in SHG 
responses seen for the relatively lightly doped samples becomes a monotonic increase. The 
same behavior is seen for both P and As doping, indicating that the effect is due to car- 
rier concentration, not type of dopand. Experiments on p-type substrates show additional 
differences, illustrating that carrier type is also important. These results will be presented 
elsewhere. 

Possible complications include bulk-quadrupole and electric-field-induced-second-harmonic 
(EFISH) contributions. [14] However, any bulk contribution will show no time dependence 
and hence, even if present, can be neglected. EFISH contributions are isotropic, result from 
charge trapping in the overlayer, and can only be a factor for heavily doped material. [15] We 
investigated EFISH for our samples by illuminating them in a glove box that was pressurized 
with N2 to prevent oxidation. For H termination the SHG intensity did not increase, as 
expected. For heavily doped oxidized samples the signal increases by about 15% before 
saturating. Thus EFISH is not a factor, and the hyperpolarizabilities that we determine are 
an accurate representation of the oxidation chemistry of the different bonds. 

To investigate whether this difference extends beyond the top two Si layers, we obtained 
thickness data on the same set of H-terminated samples as a function of time using SE. 
Details are provided elsewhere. [T6] Given the demonstrated dependence of oxidation kinetics 
on doping, the results are surprising. Except for an initial transient, Fig. 3 shows that the 
oxide thicknesses increase at a common rate up to about 1 nm, when the measurements were 
terminated. The final data in Fig. 3 include an approximately 0.2 nm contribution from 
adsorbed hydrocarbons. In fact weeks after these measurements the thicknesses remains 
almost constant. This is consistent with other data that we have obtained, including those 
of a sample that we have measured occasionally for over 30 years. These results indicate 
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that even if the initial chemical reactivities of the different types of bonds are significantly 
different, this distinction disappears after about a monolayer or so of oxidation. 

The limiting ca. 1 nm oxide thickness can be understood in terms of the electronic co- 
herence length for medium range order (MRO) in noncrystalline SiC>2. [Ej MRO extends 
beyond the first- and second-nearest-neighbor distances of ~0.16 to ~0.31 nm, respectively. 
As deduced from the position of the first sharp diffraction peak (FSDP) in the Fourier trans- 
form of X-ray diffraction data, MRO includes third- and fourth-neighbor pair correlations 
extending to 0.5 nm, as well as rings of bonded atoms. A coherence length of ~1 nm, obtained 
from the width of the FSDP, has been interpreted in terms of rigid 6-member rings encap- 
sulated by more compliant 5- and 7-member rings. [T8] This defines a characteristic cluster 
size that is effectively-strain-free. These clusters are connected through the more compliant 
lower-symmetry rings and any further degree of longer range order, e.g., nano crystalline 
periodicity, is not present. Given that the free surface provides no constraints, oxidation 
can then be expected to proceed up to a thickness where the encapsulated 6-member rings 
are formed, at which point additional energy would be required for further oxide growth. 
We thus conclude that the formation of a native oxide on Si is strain-driven to a thickness 
of ~1 nm, after which it displays long-term stability measured in years. This process is 
not dependent on the particular growth surface, as evidenced by the symmetry of the SHG 
response. 

In conclusion, our data show a distinct carrier-concentration dependence of the initial 
oxidation process of H-terminated (lll)Si. Oxidation of the outer-layer Si-Si back bonds is 
significantly enhanced with strong n-type doping. To our knowledge this is the first report 
of the dependence on carrier concentration of the chemical reactivity of different classes of 
bonds of this technologically important material. It is also the first report to recognize that 
while chemical reactivities of bonds are doping-dependent, the terminal thickness of the Si02 
overlayer is independent of carrier concentration, and is determined independently by the 
MRO length scale of Si0 2 . 



5 



[1] J. Mizsei, Vacuum, 67, 59 (2002). 

[2] H. Kakiuchi, H. Ohmi, M. Harada, H. Watanabe, and K. Yasutake, App. Phys. Lett., 90, 

091909 (2007). 
[3] H. Yamada, J. Appl. Phys, 86, 5968 (1999). 
[4] J. R. Ligenza, J. Phys. Chem., 65, 2011 (1961). 

[5] E. A. Irene and D. W. Dong, J. Electrochem. Soc, 125, 1146 (1978). 
[6] D. Fuoss and J. A. Topich, Appl. Phys. Lett., 36, 275 (1980). 

[7] J. Takizawa, S. Ohno, J. Koizumi, K. Shudo, and M. Tanaka, J. Phys.: Condens. Matter, 
18, L209 (2006). 

[8] Y. Kamiura, K. Hasegawa, Y. Sano, Y. Mizokawa, and K. Kawamoto, J. Vac. Sci. Technol., 
B 20, 2187 (2002). 

[9] W.B. Ying, Y. Mizokawa, K. Tanahashi, Y. Kamiura, M. Iida, K. Kawamoto, and W.Y. 

Yang, Thin Solid Films, 343/344, 393 (1999). 
[10] E. J. Adles and D. E. Aspnes, Phys. Rev. B, 77, 165102 (2008). 
[11] G.D. Powell, J.F. Wang, and D.E. Aspnes, Phys. Rev. B, 60, 205320 (2002). 
[12] B. Gokce, E.J. Adles, D.E. Aspnes, and K. Gundogdu, PNAS, 107, 17503 (2010). 
[13] CP. Wade and C.E.D Chidsey, Appl. Phys. Lett., 71, 1679 (1997). 
[14] H. Hirayama, K. Watanabe, and M. Kawada, Appl. Surf. Sci., 100, 460 (1996). 
[15] Z. Xu, X.F. Hu, D. Lim, J. G. Ekerdt, and M. C. Downer, J. Vac. Sci. Technol. B, 15, 1059 

(1997). 

[16] D.E. Aspnes and A. A. Studna, Appl. Opt., 14, 220 (1979). 

[17] G. Lucovsky, Phys. Stat. Sol. A., 207, 631 (2010). 

[18] J. Du and L. R. Corrales, Phys. Rev. B, 72, 09221 (2005). 



6 



I. FIGURE CAPTIONS 



Figure 1: (a) Evolution of SHGA during air exposure of a H-terminated P-doped (lll)Si 
sample with a carrier concentration n = 10 12 cm" 3 , (b) As (a) but for n = 10 14 cm" 3 . 
(c) As (a) but for n = 10 16 cm -3 , (d) Evolution of the average of these SHGA signals at 
azimuth angles of 0, 120, and 240°. Data for the As-doped sample with n = 10 18 cm" 3 are 
also included. 

Figure 2: Evolutions of the hyperpolarizabilities of the up- and back bonds of these 
samples as calculated in the ABM. 

Figure 3: Evolutions of oxide thicknesses for n-type Si as measured by SE. 
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FIG. 2. 
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